Is Cosmology Compatible with Blue Gravity Waves ? 
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A primordial gravitational wave background with positive (blue) spectral index is expected in sev- 
eral non-standard inflationary cosmologies where the stress-energy tensor violates the null energy 
condition. Here we show that a sizable amount of blue gravitational waves is compatible with 
current cosmological and astrophysical data. So far most of the works on parameter estimation 
from cosmic microwave background data have assumed a negative or negligible spectral index. The 
present limits on cosmological parameters, especially on the scalar spectral index, widen up consid- 
erably when one allows also for blue tilts of the tensor spectrum. Since the amplitude of the CMB 
B-mode polarization is larger in these models, future data from Planck are likely to provide crucial 
measurements. 
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I. INTRODUCTION 

During the last decade or so, the fluctuations and the 
polarization of the cosmic microwave background (CMB) 
have proven to be the most valuable observational tool 
to constrain cosmological models (see e.g. ||l[). To a 
big extent this is due to that fact that CMB physics is 
sufficiently simple so that it can be calculated to high ac- 
curacy with moderate computational investment. There- 
fore, accurate experimental results can be compared with 
accurate theoretical predictions. CMB observations are 
so valuable, since they provide a window to the physics 
at inflation (see e.g. [| |, |, | |, |, |, |ll|, |l|). 

This physics most probably involves the highest ener- 
gies ever probed by an experiment, more the ten orders 
of magnitude higher than energies achieved in terrestrial 
experiments like LHC at CERN. An important predic- 
tion of simple inflationary models is the production of 
a stochastic background of gravity waves (||l^) with a 
slightly tilted spectrum. 
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(1) 



where e = —H/H^ denotes a slow roll parameter from 
inflation and H is the Hubble rate during the inflation- 
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ary stage. Since in standard inflation e is strictly posi- 
tive p5| , in the usual parameter estimation routines, the 
tensor spectral index is assumed to be "red" {ut < 0, 
see e.g. |ll[) or negligible However, this does 
not rule out a priori the possibility that the spectral 
index of tensor modes might be positive, ??t > 0, i.e. 
"blue" . For instance, in the string gas cosmology set-up 
Ip^ , where scalar metric perturbations are thought to 
orig inate from initial string thermodynamic fluctuations 
||l7| , a spectrum of blue gravity waves (BGW hereafter) 
is predicted The latter is also a generic prediction 
of a class of four-dimensional models characterized by a 
bouncing phase of the universe. To induce the bounce, 
the stress-energy tensor must violate the null energy con- 
dition (NEC). In a spatially flat FRW metric, the NEC 
corresponds to the inequality H < and is ultimately 
responsible for the red tensor spectrum in standard in- 
flation. In the class of bouncing models ||l^ the scalar 
metric perturbations are originally of isocurvature nature 
and they are subsequently transformed into adiabatic 
ones. The violation of the NEC allows a BGW spec- 
trum. The same is true in the so-called super-inflation 
models |20| where inflation is driven by a component vi- 
olating the NEC. BGW may occur also in scalar-tensor 
theories and f{R) gravity theories. 

While we are aware that all the scenarios mentioned so 
far are not free from criticisms | pl| , due to our igno- 
rance of the dynamics of the very early universe, in this 
paper we will assume a more phenomenological attitude 
accepting the possibility that tensor modes might have a 
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blue spectrum. This option, together with the rather sur- 
prising fact that, as far as we know, this parameter range 
has not yet been fully analyzed, prompted us to investi- 
gate whether the presence of such a BGW spectrum is 
compatible with current cosmological and astrophysical 
observations^. Moreover, a BGW spectrum can affect 
the constraints on other cosmological parameters which 
have recently been derived from present observations of 
the cosmic microwave background. It is not surprising 
that the limit on the tensor to scalar ratio r depends on 
this assumptions, but as we shall see below, also other 
quantities like the scalar spectral index, rig, the physi- 
cal baryon density Vl^h'^ = ujh and the cold dark matter 



density, 



do. (Here ilt is the baryon density 



parameter, f2c is the cold dark matter density parame- 
ter and h is the value of the Hubble constant in units of 
lOOkm/sGc/Mpc.) 

A careful study of the current astrophysical constraints 
on the spectral index of the gravity wave background 
has been recently presented in |23). If we assume that 
the stochastic gravity wave spectrum which affects CMB 
fluctuations, at wave lengths of the order of the Hubble 
scale from Hq^ ~ 3 x 10^''ft.~^scc to about 10~^Hq^, 
extends up to wavelengths of about lO^sec relevant for 
the timing of millisecond pulsars, constraints of the or- 
der of riT ^ 0.53 can be found Similar 
constraints can be obtained from the LIGO interferom- 
eter (see |2^) while applying the nucleosynthesis bound 
on a gravity wave background yields the best constraint, 
riT ^ 0.15. However, these limits have been obtained 
by assuming the CMB upper limit on the scalar/tensor 
ratio of primordial perturbations (taken at wavelength 
k = O.Q02hMpc~^) r < 0.3 taken from the recent 
WMAP-I-SDSS analysis of I) which assumes nx = 0. 
Since a correlation clearly exists between the tensor am- 
plitude and the spectral index (e.g., clearly no constraint 
on ht can be derived if r is negligible), here we provide a 
further analysis by properly analyzing the full CMB data 
and correlations. Moreover, those limits apply only if the 
rather bold extrapolation is made that the gravity wave 
spectrum has a constant spectral slope nx over the range 
of many orders of magnitude. The wavelength relevant 
for LIGO is ~ lOsec and for nucleosynthesis even a fixed 
spectral index up to the Planck scale, ^ 10~''^sec, is as- 
sumed to obtain the above rather stringent limits ||2^, ^ . 
Since these extrapolations are so bold (we are not aware 
of any physical situation where a scaling behavior ex- 
tends over more than ten orders of magnitude), we shall 
also study the case where the limits derived in Ref. ||2^ 
do not apply. 

Below we also show how the limits on other parame- 
ters are affected when we allow for BGW spectra. We 
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FIG. 1: 68% and 95% likelihood contour plots from WMAP 
data where no external prior on the spectral index of the 
gravity wave background is assumed. 
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^ A notable exception is the paper |13[ . In this work the authors 
consiilered a BGW but with an upper Umit tit < 0.2 which, as 
we will see, is much smaller than the range of values studied here. 



FIG. 2: 95% likelihood contour plots on the riT-r plane from 
WMAP data when external prior on the spectral index of the 
gravity wave background are assumed. 



investigate the 3-year WMAP data ^jj and we produce 
forecasts for the very near future Planck satellite ex- 
periment. 



II. ANALYSIS METHOD AND RESULTS 

The method we adopt is based on the publicly avail- 
able Markov Chain Monte Carlo package cosmomc 
with a convergence diagnostics done through the Gel- 
man and Rubin statistics. We sample the following eight- 
dimensional set of cosmological parameters, adopting flat 
priors on them: the baryon and cold dark matter densi- 
ties LUb and uJc, the ratio of the sound horizon to the angu- 
lar diameter distance at decoupling, 9s, the scalar spec- 
tral index ns, the overall normalization of the spectrum 
A at k = 0.002 Mpc~^, the optical depth to rcionization, 
T, the tensor-to-scalar ratio r at fc = 0.002 Mpc~^ and. 
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FIG. 3: 95% likelihood contour plots on the ns-r plane from 
WMAP data when external prior on the spectral index of the 
gravity wave background are assumed. 
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jiT < 0.04771og( — -— — j PULSAR, (2) 
nr < 0.0223 log ( ^"^^ ^ LIGO, (3) 



riT < 0.00714 log 



3.4 X lO^h^ 



BBN. (4) 



It is important to notice the dependence of the LIGO 
and BBN constraints on h. 

Ill addition to this astrophysical constraints we have 
also simply considered a prior — 1 < ?it < 20 without 
taking the above limits into account, since it may very 
well be that the gravity wave spectrum is modified on 
some scale below those needed for the limits of Eqs. (H) 
to (^. The results are also compared with a 2 — cr upper 
limit of JiT < 

A first analysis is made in the theoretical framework 
described above with no upper limits on ut- The likeli- 
hood contour plot in the riT-r plane is reported in Fig- 
ure ^ When only the WMAP data is considered, no 
upper limit on ut alone can be derived. The reason is 
simple: gravity waves with BGW spectra can always be 
accommodated with the WMAP data by lowering the 
tensor/scalar ratio r. However it is possible to derive 
the following 95% upper limits on ut in function of r: 
riT < 2.1 for r > 0.1 and nr < 3.2 for r > 0.01. In the 
range 0.05 < r < 0.6 and 0.75 < nr < 3 we derive the 
following upper limit fit at 95% c.l. from WMAP only: 



r + 0.68nT-0.l2n^<l (5) 
while including the ACBAR data we obtain: 



FIG. 4: 95% likelihood contour plots on the ujb-coc plane 
from WMAP data when external priors on tit are assumed. 
The sohd line corresponds to the WMAP-I-PULSAR, the 
dashed hue to WMAP-I-LIGO, the long-short dashed to 
WMAP-I-BBN and the dotted line to WMAP+nr < case. 



finally, the tensor spectral index tit- Furthermore, we 
consider purely adiabatic initial conditions and we im- 
pose spatial flatness. 

We include the three-year WMAP data ||l| (tempera- 
ture and polarization) with the routine for computing the 
likelihood supplied by the WMAP team. Moreover we 
also include in a separate analysis the new, high quality, 
fine-scale measurements from the ACBAR experiment 

(0). 

We have studied three different possible priors on the 
combination between ut and r as derived in Rcf. 
from pulsar timing, LIGO and nucleosynthesis with the 
Planck scale as the upper cutoff: 



r -f 0.41nT - 0.06n^ < 0.7 (6) 

again at 95% c.l.. Assuming r > 0.05 and no prior on 
nr, we found from WMAP+ACBAR -1.51 < tit < 2.62 
at 95%c.L 

These limits, while new and on very different scales, 
are however rather weak if compared with the upper lim- 
its described in Eqs. ^ to (||) above. It is however in- 
teresting to notice that even when these stronger limits 
are considered, BGW spectra are still compatible with 
CMB observations. Since this contribution has not been 
usually considered in CMB parameter analysis, it is inter- 
esting to study the possible effect of such a BGW back- 
ground. 

In Figure ^ |3| and ^ we show how cosmological pa- 
rameters as inferred from the WMAP 3-ycar data are 
affected by a BGW spectrum. 

In Figure ^ we see how current external priors on ut 
affect the WMAP bounds on r. As we can see, current 
bounds on r obtained for ut < can be relaxed by as 
much as a factor 1.4 when BGW spectra in agreement 
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with LIGO and PULSAR data are assumed. In partic- 
ular we found the following 95% c.l. upper limits on 
r: r < 0.64, r < 0.65, r < 0.55 and r < 0.47 when 
the WMAP data is combined with the PULSAR, LIGO, 
BBN and < priors respectively. We found that the 
constraints with no external prior such that ut < 1 are 
very similar to the PULSAR and LIGO cases. 

In Figure ^ we see how the parameter ranges for the 
pair {ns, r) are affected by a blue tilt in the tensor spec- 
trum. First, it is interesting to note that a larger tensor 
contribution is allowed and the scalar spectral index can 
then be bluer. This comes from the well known par- 
tial degeneracy between rising the scalar spectral index 
which reduces the Sachs Wolfe plateau if the height of the 
first acoustic peak is fixed, and at the same time rising 
the tensor contribution which enhances the Sachs Wolfe 
plateau. This degeneracy is strengthened if > since 
then tensor fluctuations in the CMB extend to somewhat 
smaller scales. Disregarding the very speculative nucle- 
osynthesis prior, the 2a limit on the scalar spectral index 
is enhanced from = 0.981 ±0.026 to = 0.994±0.032 
when allowing for BGW spectra consistent with pulsar 
timing. BGW's affect also the constraints on the optical 
depth parameter, which moves from r = 0.091 ±0.030 to 
T = 0.097 ± 0.032 for the WMAP+PULSAR case. 

As one sees in Figure |[ enhancing the scalar spectral 
index is slightly correlated with enhancing Wf, which leads 
to stronger Silk damping. The higher baryon density is 
compensated by a lower cold dark matter density. This 
not in a way so that uim = iVc+tOb would remain constant, 
but so that the total matter density decreases which also 
reduces the acoustic peaks, to compensate for the blue 
spectral index. The 2a limits on cUc and ut change from 
LUb = 0.0228 ± 0.0009 and lUc = 0.100 ± 0.009 to ujb = 
0.0232 ± 0.0011 and Uc = 0.098 ± 0.01 respectively when 
allowing for blue tensor spectra consistent with pulsar 
timing. The effect of BGW's is therefore smaller on those 
parameters but still measurable. 

While current CMB constraints are clearly affected by 
the inclusion of a blue gravity wave component it is inter- 
esting to investigate the impact on future measurements 
as those expected from the Planck satellite. The first 
important observation is that a blue tensor component 
dramatically enhances the probability for the detection 
of a gravity wave background with PLANCK. In Fig- 
ure H we plot the best fit angular anisotropy power spec- 
tra to the WMAP data for the different external priors 
on ut- As we can see, all the models give nearly iden- 
tical power spectra, i.e. it is impossible to discriminate 
between those model with current CMB observations. 

CMB polarization is, however, a very promising tool 
for detecting a BGW background. The statistical prop- 
erties of CMB linear polarization are fully characterized 
by two sets of spin-2 multipole moments with opposite 
parities |30|. As well known in the literature (see e.g. 
pl| , ^ ) the magnetic- type modes (B or curl modes) are 
produced by tensor metric perturbations and not scalar 
perturbations. A detection of B-mode polarization would 
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FIG. 5: Top Panel: CMB Temperature anisotropy angular 
power spectra for the best fit models to current WMAP data 
for diflFerent choices of external priors on nr. Bottom Panel: 
Contribution to the B mode polarization. 



thus provide good evidence for a primordial background 
of gravity waves. 

As we can see from the bottom panel of Figure I the 
expected amplitude of the polarization B mode differs 
dramatically between the models. In particular, it is im- 
portant to notice that a blue spectral index could en- 
hance the B mode power at ^ ^ 100 by several orders of 
magnitude if compared with models with ut < 0. Also 
in the plot we also report the expected sensitivity from 
the Planck experiment (dashed horizontal line). Most 
of the BGW spectra models are above the sensitivity of 
PLANCK. 

We have simulated future eosmological data with a 
noise spectrum based on the expected Planck configura- 
tion. In particular we assume an experimental sensitivity 

— 1/2 

Wp = 81/iif, beam size Ofwhm ~ 7.1' and sky cov- 
erage fsky = 0.8 We also assume a spatially flat model 
with parameters h = 0.7, ojc = 0.120, LUb = 0.022, ns = 1 
and r = 0. We find that the constraints on tensor modes 
expected for this model will be equally affected. The 2a 
limit on r is indeed relaxed by ~ 30% from r < 0.025 to 
r < 0.035 when allowing for blue tensor spectra consis- 
tent with pulsar timing. 



III. CONCLUSIONS 

In this paper we have analyzed to which extent a pri- 
mordial background of gravitational waves with positive 
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spectral index is compatible with current observations. 
We have found that a considerable part of the param- 
eter space of cosmological models is in agreement with 
current WMAP data and with upper limits coming from 
the LIGO experiment and from pulsar timing. Several 
non standard models can produce such a background and 
there is therefore no strong theoretical reason to exclude 
these models in current parameter estimation analysis. 
With this in mind, we have shown that, if the nx < 
assumption is relaxed, also parameters which are not di- 
rectly related to the tensor component are affected. Es- 
pecially, the scalar spectral index can be significantly 
bluer than in models which allow only for a red tilt of 
the tensor spectrum. Moreover, we have found a smaller 
but significant effect on the baryon and cold dark matter 
density constraints. This example shows once more that 
when assuming any limits on cosmological parameters it 
is very important to be aware of the model assumptions 
which went into their derivation. 

Future Planck estimates can also be affected by BGW 



spectra, especially when constraints on r are considered. 
However they will be nearly entirely free of the additional 
degeneracies which are still seen in the WMAP data. As 
a final remark, we like to point out that blue spectra are 
able to produce significantly larger contribution to the 
_B-modc polarization spectrum. An excess of B-mode 
polarization at £ ~ f 00 would therefore provide convinc- 
ing evidence for the BGW models investigated here. 
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